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I. INTRODUCTION 



For a number of years black-hole (BH) quasi-normal modes (QNMs) have been of great interest to both gravi- 
tational theorists and gravitational wave experimentalists. QNMs have been studied in the contexts of ringing of 
astrophysical BHs where QNMs are relevant to gravitational wave astronomy [IHS], loop gravity interpretations [3], 
lower dimensional BHs [5], and in the AdS/CFT correspondence [6] where BHs are thermal states in the field theory 
and their decays correspond to a decay of a perturbation of this state [7]. As such, knowledge of QNMs will provide an 
understanding of the stability of BHs against external perturbations; where more recently the interest in the Hawk- 
ing radiation and QNMs of BHs has stemmed from the idea that rmm-TeV scale BHs might be created in particle 
accelerators such as the Large Hadron Collider (LHC) ^8-12^. QNM oscillations have been found in perturbation 
calculations of particles falling into Schwarzschild [T3] and Kerr BHs [T3] and in the collapse of a star to form a black 
hole [H]. 

Let us first recall that QNMs are the damped, resonant, non-radial perturbations of BHs which are characterized 
by a spectrum of discrete, complex frequencies, with real part determining the oscillation frequency, and an imaginary 
part that determines the rate at which each mode is damped as a result of an emission. For a given kind of physical 
perturbation the complex frequencies are uniquely determined by the mass and angular momentum of the BH, the 
angular harmonic indices (/, m) of the deformation, and the degree of the harmonic of the mode n. As has been shown 
in many contexts, the fundamental equations describing the perturbations of BHs can be reduced to a single second- 
order ordinary differential equation similar to the one-dimensional Schrodinger equation for a particle encountering 
a potential barrier. However, the nature of the potential precludes an exact, closed-form solution in terms of known 
functions. That is, our equations for wave propagation reduce to 

^ + (.^-n.))V = o, (1.1) 

where x is a spatial co-ordinate. The BH event-horizon is mapped to — oo in the radial co-ordinate, and the potential 
V{x) is a function of the spatial co-ordinate only. Note that we have proceeded in the usual way by assuming a time 
dependence, '^{t,x) = e*"'?A(a;). 



This separation is not restrictive since once we have a solution for Eq.(l.l), a genera l tim e dependent solution can 
be given as a continuous Fourier transform of such a solution, where the form of Eq.(l.l) parallels a normal mode 
analysis. We shall now restrict ourselves to asymptotically flat space-times such that the potential V{x) is positive 
and satisfies 

V{x) ^ , X'^±oo . (1.2) 
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Note that there exist no normal mode expansions for such systems since such potentials do not allow bound states. 



From the form of the potential in Eq.(1.2| we have that near to the boundaries ±00 the solutions behave as plane 
waves, 

tfj{x) ~ e=^'"^, X ^ ±00 . (1.3) 
Here in-going at the horizon means entering into the BH, whilst the QNM frequencies and the associated wave- 



functions ip (which are the solutions to Eq.(l.l)) are the QNMs specified by the boundary conditions of Eq.(1.3). 

The study of BH perturbations was initiated by Regge and Wheeler \W and extended by Zcrilli [17 . Following the 
work by Vishveshwara Chandrasekhar and Detweiler [15] pioneered numerical methods for the study of QNMs. 
More recently, semi-analytic methods have proven useful in the study of BH perturbations with the WKB(J) method 
of Schutz, Iyer and Will [201 '21], and more recently the Asymptotic Iteration Method (AIM) of Ciftci, Hall and Saad 
[22j . which is a valuable tool for computing solutions to general second order ordinary differential equations. 

In this note we extend to sixth order in the semi-analytic approximation the works of Seidel and Iyer [33] , Konoplya 
[24| . Iyer and Will [5T], and following the works of Refs. we extend to the fermion sector the Kerr BH 

QNMs. Here we shall compare the results of direct numerical, WKB(J) approximation and the AIM computation of 
QNMs of fermionic BH excitations in the Kerr space-time. Note that we shall assign c = h = G ^ 1 and M — 1 
for the calculated QNM frequencies in this work, where this note is organised as follows: In Sec[n] we shall give a 
brief overview of the Newman-Penrose (NP) formalism [27] and its application to the QNMs of massless spin-^ fields. 
Sec|III| and |IV| contain reviews of the semi-analytic methods employed here, whilst SecjV] contains the numerical 



results of the QNMs for spin-i field in a Kerr space-time, followed by a short discussion of these results. 



II. NEWMAN PENROSE FORMALISM 



In order to develop a master equation as discussed in the previous section, for a fermion, we shall introduce the NP 
formalism |27J. Let {Ai,g) be a pseudo-Riemannian manifold with metric g of dimension n, and denote the tangent 
space at p € by TpM. A natural basis for TpM is the set of co-ordinate partial derivatives = d^, although we 
may choose a basis as we wish. Therefore, we choose one which is orthonormal with the signature of the metric g on 
M, that is 

ff(ea,e6) = riab , 

where r]ab is a block diagonal matrix in the NP formalism. A vector may be expressed as a linear combination of basis 
vectors such that 

e,, = e^e, , (2.1) 
and the operators are n x n invertable matrices which satisfy, 

e'-.e: = 5; , 

where Latin letters index tetrad co-ordinates and Greek letters index manifold co-ordinates. The operator ej^, now 
called a vielhein, is a projection operator from the co-ordinate basis to the non-co-ordinate basis. The NP formalism 
is a special case for the vielhein^ or tetrad, formalism in which tensors defined on a space-time are projected on to a 
complete basis vector set at each point p. This basis set is the set of complex null vectors 1, n, m and m, of which 1 
and n are real and m and m are complex conjugate. By the null condition of this basis vector set, it is clear that 

ll = nn==mm = rnm = 0. 



In addition, these basis vectors satisfy the orthogonality relations, 

lm = lm = nm = nm = 0, 
as well as the normalisation conditions 1 • n = 1 and m • m = —1. The corresponding fundamental matrix rjab is 






1 








1 




















-1 








-1 
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to which there corresponds a basis, ei = 1, e2 = n, 63 = m and 64 = m. When used to define directional derivatives 
we assign to each basis vector a special symbol, ei = D, 62 = D*, 63 = 6, and 64 = S* . 

Having developed the NP formalism, we begin our analysis by considering the Dirac equation in the null tetrad 
basis, 



{-ijf'S/f, + to) = 



where is the covariant derivative taking the form 



The spin connection 0;^"^ satisfies the vielbein ansatz, 

da I a b j^a a r\ 



(2.2) 
(2.3) 
(2.4) 



and 7" — — with 7" being the a Dirac matrix. Note that the spin connection w^"^ may be rewritten in 
terms of the affine connection as 



The spin connection becomes the affine metric connection under the following co-ordinate transformation. 



We shall focus our attention on the massless limit of Eq.(2.2), 

7^V^* = , 



rewriting ^I' in terms of Weyl spinnors and ■(/;, that is, 

^ = 



and substituting Eq.(2.8l into Eq.(2.7) yields the dynamical Weyl spinor equations, 

'^^bV^^'^^O, and a^.V^^'' = 0, 



(2.5) 
(2.6) 
(2.7) 

(2.8) 
(2.9) 



in the chiral representation, cr" and ct" are Pauli matrices, and and V/i are the the covariant derivatives of ■0° and 
T/;", respectively. The Pauli matrices take on new forms in the vielbein formalism for curved space-time, p'' = e^cr'^, 
= ej^fja, = ^a^'^ P^^ = ^li^a, where are the vielbein basis elements and cr^ are Pauli matrices. Explicitly 
we find 



P ah 



and 



ah 



Following the approach of Ref. [28], the covariant derivative of Eqs.(2.7) and (2.9) may be written as 

P\i>'^f^r = 0, and p^V^V^^ = 0. 
Consider now the 4-dimensional Kerr background space-time with line element 



(2.10) 



(2.11) 



ds^ = - 1 - 



2Afr\ 



AMra sin^ 9 
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dtde 



, ? / 9\ 17 ? 2Mrasin 

dr^ + dO^ + ( + + — 



sm 



(2.12) 



where p^ ^ r"^ + cos^ 6, A 
define p — r + ia cos 9 and p* 
geometry as 



+ — 2Mr and a = -p^ e [0; 1) is the rotation parameter of the BH. It is useful to 
r — ia cos 9 such that p^ — pp* . The basis of the NP formalism as defined on the Kerr 



Ij = — (A, 0, -aA sin^ ( 



V2p 



— lasmf 



sm( 



1 

rii — — - 
' 2p' 

•ffiA — to! 



(A,p^O, -aA sin^ 
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and 



\ f i 
iasinO, 0, 1, 



K 


= '^Sll , 


P 


= 1^314 , 


e 


a 


= "^313 , 


A* 


= ^243 , 


7 


A 


— ^^244 J 


T 


— "^312 , 


a 


V 


— ^242 , 


TT 


= 1^241 , 





\/2p V ' ' ' sin ( 

We may now define the Ricci rotation coefficients (spin coefficients in that which foUows) with respect to the curved 
space-time, 

^ (W211 + 1^34l) , 

^(c.212+^324) , (2.13) 

1 (.^213 + "^344) , 

2 ('^213 + 1^343) • 

Recall that the complex conjugate of any spin coefficient is obtained by replacing the numerical index 3 with the value 
4 and vice versa [28]. The resulting equations are 

(£> + e - p) V'" + ((5* + TT - a) = , 
(£>* + M - 7) V^^ + (5 + /3 - r) V" = 0, 

{D + e* -p*)i!^ + {6 + TT* ~a*)i^ = 0, 

{D* + ^i* + {6* + P* -T*)i!^ = 0. 

It is clear that basis vectors are mapped to derivative operators when applied to the tangent space as tangent 
vectors to the wave-functions ip{t,r,9,(l)) — e*'^"*+™'^''(/'(r, 0). Therefore, we write 

l = D = Vo , n^D* ^-^Vl, 

m = S = -^Cl , m = (5* = -Co , 

V2p V2p* 



where 



^ n K ^ r-M . ^ K ^ r-M 

n ina sin 6 ^+ n ^ " sin ( 

Cn^ Oe + Q + 7: z , Cl^ = O0 -Q 



tand p 

with K = (r^ + a^) uj -I- am, and Q = aw sin H -. 

^ ' sm 6* 

The spin coefficients are linear combinations of basis vectors in the NP formalism. We may therefore express 

the collection of spin connections in terms metric tensor elements. We shall choose the representation of the spin 

connections so as to simplify the algebra. A suitable choice for this assignment of elements is 

1 ^ iasiaO 

^ ~ p* ' 2\/2p*tane ' ^ ~ \/2{n*\'^ 



p- 

iasinfl ^ r — M 

a — TT — P* , 



which results in. 



K = (T = X = iy = e = 0. 

Now let /i = /0*V'°, /2 — V'^: 91 = V"^; .92 = —pip'^- Followiug these assignments, it can be shown P5| that the Dirac 
equations reduce to 

i^oA + ^/:,/2 



AV\f2-V2C\h 

2 2 



0, 


(2.14a) 


0, 


(2.14b) 


0, 


(2.14c) 


. 


(2.14d) 
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Using Eq.(2.14|, we separate the Dirac equation into radial and angular parts by making the ansatz: 



Substituting Eq.(2.15l into Eq.(2.14) yields, 



and 



CiSi = — AS'_i , 

2 2 2 

= xsi , 

2 2 2 



(2.15a) 
(2.15b) 
(2.15c) 
(2.15d) 



(2.16a) 
(2.16b) 

(2.17a) 
(2.17b) 



where R' i = %/2i?i . 



Decoupling Eq.(2.16l yields, 



5^1 = , 



^2 ^2 

where Si is the adjoint of S'_i under the replacement oi 9 hy t: — 9. Similarly, decoupling Eq.(2.17) yields. 

Apt Pol = , 



(2.18) 



(2.19) 



and V ARi satisfies the complex conjugate relation. Once decoupled, Eqs.(2.18) and (2.19) are eigenvalue equations 
with known eigenvalues A and unknown eigen-functions R^. i and i . 

For the purpose of computing QNMs we shall focus our interest on the radial equation for the spin-^ field. We map 
Eq.(2.19l to the tortoise co-ordinate x, 



A 

K 



(2.20) 



K 



where K — — eliminates the co-ordinate singularity at the horizon. After this co-ordinate transformation, operators 



I?o and I?o simplify to, 



K 



A \dx 



, t K f d 
anaVk — -r [ -, 



(2.21a) 
(2.21b) 



The assignments Pi = V~ARi and P_i = in conjunction with Eq.(2.21), lead to the new form of the radial 

equations 



dx 



iuj Pi 



dx 



iuj\P_i 



Va 

A^P. 
K 5 



(2.22a) 
(2.22b) 



If we define, 



Pi ±P 



(2.23) 
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and combine Eqs.(2.22a) and (2.22b) we obtain 



which may be separated to yield 



with 




(2.24) 



^ ' dx \ K 



(2.25) 



Eq.(2.24) is now in the form of Eq.( |l.l[ ) with a potential given by Eq.(2.25). Eq.(2.24) is the QNM master equation 
for a spin-i field in a Kerr space-time. In the following two sections we shall review two semi-analytic approaches to 
solving these equations for an ultimate comparison with the direct numerical results of Refs. I30j . 



III. WKB(J) METHOD 

We shall begin with the WKB(J) approximation of Schutz, Iyer and Will [2ni [H] which is a general method for 
generating approximate solutions to ordinary linear second order differential equations and has been used extensively in 
various BH cases (see Ref. [21] and references therein). For further discussion of this method, see Refs. pnE51IMll31j . 
We may use the WKB(J) method to generate approximate solutions to differential equations of the form 

^ g(a;)* = 0, (3.1) 



dx"^ 

for a smooth function, at 0(1), Q{x) = cj^ — f{x) with oj E C Now assume that has a mode expansion of the form, 

*(x, e) = Ae^^=o ^SAx)e-^ ^ ^3 2) 



where A £ C and e <C 1 tracks the order of the expansion. By substituting Eq.(3.2) into Eq.(3.1|, and equating 
powers of e, yields 



So{x) =±ij dt ^/Qitj , (3.3) 

and 

5i(x) = -ilnQ(x), (3.4) 



where the two choices in sign in Eq.(3.3l correspond to either incoming or outgoing waves at ±cxd. When x — -|-c», 
Q{x) — > o;^ such that Sq +iujx for the outgoing wave to infinity and 5*0 — )■ —iojx for the incoming wave from 
infinity. Similarly, for x — >■ —cx), 5o — > +iujx for the wave incoming from a; — ^ — oo, while Sq — > —ujx to a wave 
outgoing to a: — ^ — oo, see Figjl] Designating these four solutions. 



III ^-\-ujx Willi „—ujx ^ V 

(3.5) 



The general solutions to regions I and III are 

* - Z'j^-^Y + Zf^^^^^ , region I , 



IN it,- , 70C/T,T,/ ^^„;„„ TTT 

(3.6) 



* - Z\fj^Jj + Zyij^ , region III 
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FIG. 1: The reduced potential barrier as used in the WKB(J) method. 



where Zj , Z/, Zjjj, Z^jj e C. 

The ampHtudes at +00 are related to amphtudes at —00 through the hnear relation 



Zjjj \ ^ / Sii Si2\ f Z[ 
Zjjj J \ S21 S22 J \ Z_ 



(3.7) 



We match the solutions in regions I and III with a solution in region II through the two turning points of Q{x) = 0. If 
the turning points are closely spaced {—Q{x)max ^ then the solution in region II can be well approximated 

by a Taylor series expansion, 

Q(x) = Q{xo) + ^Q"{xo) {x - xof + O ((x - .To)') , (3.8) 

where xq is the location of the maximum of the function Q{x) and the prime notation signifies differentiation with 
respect to x. We shall assign Qo — Q (xq) in what follows. Region II corresponds to 



We introduce new functions 



\x-xo\<J^^e'^ . (3.9) 



t = (4A:)JeT {x - xq) , (3.10) 
and.+ i = ^. (3.11) 



2 



Rewriting Eq.( 3.1 ) as 



d^fX / 1 1 



which has, is terms of parabolic cylinder functions F^{t), the general solution 

m{t) = AF^ + BF^^idt) . (3.13) 
It can be shown that the large |t| asymptotics of this solution give an S'-matrix, 



where T is the Gamma function and i? is a general exponential function depending on the nature of the potential 



Q. When expanding to higher orders, the S matrix has the same general form as shown in Eq.(3.14), with modified 
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expressions for the off-diagonal elements. Note that for a BH no waves are reflected by the horizon, such that = 0, 
and due to the QNM boundary conditions there are no waves incoming from infinity, ZL = 0. These conditions are 
satisfied in Eq.(3.7) only if r(— z^) — 0. Consequently, G Z, and 



n = 0, 1, 3, . 



(3.15) 



The above equation gives the quantisation condition which labels QNMs by an overtone number n = —v. Since the 
coefficients in the ^-matrix that connect amplitudes near the horizon with those at infinity depend only on a (related 
to the overtone number n), we need only expand Eq.(3.2 1 to higher orders in e to obtain the WKB formula to higher 
orders. The result has the form 



Ao - A.s - A4 - Ar 



Afi 



n + 



1 



(3.16) 



where the correction terms A4, A5, Ag can be found in Ref. [3T] and the included URL. 



IV. ASYMPTOTIC ITERATION METHOD 



Our second semi-analytic technique is the AIM, which produces an exact solution to an ordinary second order 
partial differential equation under the assumption that the coefficients in the Taylor expansion of the differential 
equation form an asymptotic series. To see this, consider the homogeneous linear second-order differential equation 
for the function x{^)-i 

x" = Ao(a:)x' + so(a;)x, , (4.1) 

where Ao(a;) and sq{x) are functions in C'ooio-,b). In order to find a general solution to this equation we rely on the 
symmetric structure of the right-hand side of Eq.(|4.1| [22 . If we differentiate Eq.(4.1 ) with respect to x, we find that 



x" = M{x)x + si{x)x , 
where Ai = Aq + sq + (Aq)^ and si ~ s'q + sqXq- Taking the second derivative of Eq.(4.1 1 yields 

X"" = X2{x)x' + S2{x}x , 

where A2 = A'^ + si + AqAi, and Si = s'q + sqXq. The {n + 1)*'' and the (n + 2)*'* derivatives, n — 1,2, . . are found 
iteratively, 

x("+i)=A„_i(a;)x'-Hs„_i(x)x. (4.2) 

The crucial observation in the AIM is that differentiating the above equation n times with respect to x, leaves a 
symmetric form for the right hand side: 

x("+2)=A„(x)x' + .s„(a:)x, (4.3) 

Xnix) = X'n-iix) + Sn-lix) + Ao (x) A„_ i (a;) , (4.4a) 

Sn{x) = s'„_i(x) -I- SQ{x)Xn-i{x) . (4.4b) 
S„{x) Sn^i{x) 



where. 



and 



For sufficiently large n 



Xn{x) A„_i(x) 

where from the ratio of the (n + 1)*'' and the (n + 2)*'* derivatives, we have 



m , (4.5) 



d 



y("+2) A„(x'+f^x) d . ,^ A 



lnx("+^) = Vl) = 7 ^^ = f lnx("+^' = f^. (4.6) 
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This then yields, 



X^^+^H^;) - Ci exp 



dx' 



CiA„_i exp / dx'{j3 + Aq) I , 



(4.7) 



where Ci is the integration constant and the right-hand side of Eq.(4.4) and the definition of j3{x) have been used. 

(4.8) 



Substituting this into Eq.(4.2), we obtain the first-order differential equation 

x' + /3x = Ciexpf / dx'{P + xS 



The general solution can be determined as 



X{x) = exp 



dx'P{x') 



C2 + C1 / dx' exp / dx" [Xoix") + 2/3 {x") 



(4.9) 



The integration constants, Ci and C2, can be determined by an appropriate choice of normalisation. Note that for 
the generation of exact solutions Ci — 0. Eq.(4.5) gives the sequence termination condition. 



(4.10) 



which is also the QNM "quantisation condition" as given in Ref. [32] • Note that the sequence termination condition 
is true only at a particular value n since there is no guarantee that the sequence A„ is truly asymptotic. 

Improvements to the efficiency of the AIM have been made which eliminates the need to take derivatives at each 
step in the iteration process, thereby reducing the computational overhead while also improving the overall accuracy 
of the method. This is done by reducing the AIM to a set of recursion relations which no longer involve differentiation. 
For more details on this improved AIM see Ref. [25j. 



V. RESULTS AND CONCLUDING REMARKS 



Tables |Tj |TTj |III[ |TV] and |V] constitute the main results of this paper, and contain the numerically computed QNM 
frequencies for varying values of the BH rotation parameter a = jj, along with the semi-analytic WKB(J) (to third 
and sixth order) and AIM results. In each case, and where numerical results were available, they were compared with 
the semi-analytic solutions, these numerical results were taken from Refs. [211 1301 • At this point we re- iterate that 
we have set c — G = 1 and M = 1. Note that in the absence of numerical results, QNM frequencies calculated using 
the sixth order WKB(J) method were compared to the values from the AIM. 

Unlike earlier cases, such as the case of scalar perturbations, previously published numerical values exist for the 
QNM frequencies that result from spin- half perturbations of a Kerr BH [5S1 [3D] , at least for the angular numbers of 
I = and / = 1. In comparing these numerical results, the accuracy of the WKB(J) method was worst (relative to 
the pure numerical evaluation) for parameter values {a,l,n) = (0,0,0), with approximate errors of (3.55%, 3.20%) at 
third order and (0.16%, 2.16%) at sixth order, and showed steady improvement in accuracy with increasing I and n, to 
(0.13%, 0.50%) at third order and (0.02%, < 0.01%) at sixth order. However, the AIM method experienced degraded 
accuracy with increasing I and n at a = 0, ranging from (< 0.01%, 0.10%) at ^ = n = to (1.46%, 4.32%) at I — n = 2. 
However, the accuracy of the AIM was at its highest when a was at its lowest, n was at its lowest and / was at its 
highest. Note that this for an iterative of only 15, though with this iterative depth the AIM was still significantly 
slower that the sixth order WKB(J) method. 

For the QNM frequencies when I = 2 and I — 3, the AIM values were then compared to the sixth order WKB(J) 
method values. In these cases the accuracy of the AIM, as compared to the sixth order WKB(J) method was at its 
highest when n was at its lowest and the I values were higher, such as (a, I, n) = (0, 3, 0), (0.20, 2, 0), (0.40, 2 or 3, 0), 
(0.60,2 or 3,0) and (0.80,3,0). Once again, numerical processors did evaluate significantly faster for the sixth order 
WKB(J) method than for the AIM. 

In conclusion, we have, using the NP formalism, developed the Dirac equation in a Kerr background in a peda- 
gogical manner. Using this result we have applied two semi-analytic techniques to calculate the QNM frequencies for 
comparison with available numerical results [211 [3D] finding that the computational overhead for the WKB(J) method 
is lower at low n than for the AIM. As such the WKB(J) method executes considerably more quickly than the AIM. 
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TABLE I: The QNM Frequencies for the Spin-^ Perturbation of the Kerr BH with rotation parameter a = 0.00, whieh is the 
Schwarzschild limit. 



1 


n 


Numerical 


Ihird Order WKB 


Sixth Order WKB 


AIM 








0.1830 - 0.09701 


0.1765 - O.lOOli 


0.1827 - 0.0949i 


0.1830 - 0.09691 








(-6.00/0 , o.ZU/o) 


(-O.ibyo , -2.ib7oj 


(<O.Oi7o , -O.iOTo) 


1 





0.3800 - 0.0964i 


0.3786 - 0.0965i 


0.3801 - 0.0964i 


0.3800 - 0.0964i 








(-0.37% , 0.107o) 


(0.03% , <0.01%)) 


(<0.0l7o , <0.01%)) 




1 


0.3558 - 0.2975i 


0.3536 - 0.2987i 


0.3559 - 0.2973i 


0.3568 - 0.2976i 








(-0.62% , 0.40%) 


(0.03% , -0.07%) 


(0.28% , 0.03%)) 


2 





0.5741 - 0.0963i 


0.5737 - 0.0963i 


0.5741 - 0.0963i 


0.5741 - 0.0963i 








(-0.07%) , <0.01%>) 


(<0.0l7o , <0.0l7o) 


(<0.0l7o , <0.01%)) 




1 


0.5570 - 0.2927i 


0.5562 - 0.2930i 


0.5570 - 0.2927i 


0.5573 - 0.2928i 








(-0.14% , 0.10%)) 


(<0.0l7o , <0.0l7o) 


(0.05%) , 0.03%)) 




2 


0.5266 - 0.4997i 


0.5273 - 0.4972i 


0.5265 - 0.4997i 


0.5189 - 0.5213i 








(0.1370 , -0.50%) 


(-0.02%) , <0.0l7o) 


(-1.46%) , 4.32%)) 


3 







0.7672 - 0.0963i 


0.7674 - 0.0963i 


0.7674 - 0.0963i 
















1 




0.7540 - 0.2910i 


0.7543 - 0.2910i 


0.7544 - 0.2910i 












(0.01% , <0.01%) 




2 




0.7305 - 0.4909i 


0.7298 - 0.4919i 


0.7267 - 0.4928i 












(-0.42% , 0.18%) 




3 




0.6999 - 0.6957i 


0.6967 - 0.7023i 
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TABLE II: The QNM Frequencies for the Spin-^ Perturbation of the Kerr BH with a rotation parameter a — 0.20. 



1 


n 


Numerical 


Third Order WKB 


Sixth Order WKB 


AIM 








0.1836 - 0.0967i 


0.1773 - 0.0997i 


0.1837 - 0.09421 


0.1835 - 0.09651 
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0.3796 - 0.0962i 


0.3810 - 0.09601 


0.3809 - 0.09601 








(-0.39% , 0.10%) 


(-0.03% , -0.10%) 


(-0.05% , -0.10%) 




1 


0.3572 - 0.2964i 


0.3550 - 0.2974i 


0.3572 - 0.29691 


0.3583 - 0.29621 








(-0.62% , 0.34%) 


(<0.01% , 0.17%) 


(0.31% , -0.07%) 


2 







0.5750 - 0.0960i 


0.5754 - 0.09591 


0.5754 - 0.09591 












(<0.01% , <0.01%) 




1 




0.5578 - 0.2918i 


0.5586 - 0.29151 


0.5589 - 0.29161 












(0.05% , 0.03%) 




2 




0.5294 - 0. 49511 


0.5287 - 0.49751 




3 







0.7689 - 0.09591 


0.7691 - 0.09501 


0.7691 - 0.09591 












(<0.01% , 0.95%) 
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TABLE III: The QNM Frequencies for the Spin-^ Perturbation of the Kerr BH with a rotation parameter a = 0.40. 



1 


n 


Numerical 


Third Order WKB 


Sixth Order WKB 


AIM 








0.1854 - 0.0956i 


0.1798 - 0.0982i 

(-3.02% , 2.72%) 


0.1844 - 0.09321 

(-0.54% , -2.51%) 


0.1851 - 0.09511 

(-0.16% , -0.52%) 


1 



1 


0.3843 - 0.0951i 
0.3614 - 0.2930i 


0.3825 - 0.09491 
(-0.47% , -0.21%) 
0.3592 - 0.29311 
(-0.61% , 0.03%) 


0.3838 - 0.09481 
(-0.13% , -0.32%) 

0.3612 - 0.29181 
(-0.06% , -0.41%) 


0.3837 - 0.09481 
(-0.16% , -0.32%) 
0.3630 - 0.29161 
(0.44% , -0.48%) 


2 



1 
2 




0.5791 - 0.09481 
0.5628 - 0.2880i 
0.5357 - 0.48831 


0.5795 - 0.09471 
0.5636 - 0.28781 
0.5352 - 0.49051 


0.5795 - 0.09471 
(<0.01% , <0.01%) 
0.5640 - 0.28771 
(0.07% , 0.03%) 
0.5472 - 0.52041 
(2.24% , 6.10%) 


3 



1 
2 
3 




0.7744 - 0.09471 
0.7620 - 0.28621 
0.7400 - 0.48251 
0.7114 - 0.68351 


0.7745 - 0.09471 
0.7624 - 0.28621 
0.7395 - 0.48341 
0.7086 - 0.68931 


0.7745 - 0.09471 
(<0.01% , <0.01%) 
0.7625 - 0.28621 
(0.01% , <0.01%) 
0.7358 - 0.48581 
(-0.50% , 0.50%) 



TABLE IV: The QNM Frequencies for the Spin-^ Perturbation of the Kerr BH with a rotation parameter a = 0.60. 



1 


n 


Numerical 


Third Order WKB 


Sixth Order WKB 


AIM 








0.1885 - 0.09341 


0.1838 - 0.09531 


0.1870 - 0.09091 


0.1879 - 0.09251 








(-2.49% , 2.03%) 


(-0.80% , -2.68%) 


(-0.32% , -0.96%) 


1 





0.3901 - 0.09311 


0.3877 - 0.09251 


0.3888 - 0.09241 


0.3888 - 0.09241 








(-0.62% , -0.64%) 


(-0.33% , -0.75%) 


(-0.33% , -0.75%) 




1 


0.3687 - 0.28611 


0.3662 - 0.28491 


0.3682 - 0.28371 


0.3714 - 0.28111 








(-0.68% , -0.42%) 


(-0.13% , -0.84%) 


(0.73% , -1.75%) 


2 







0.5866 - 0.09241 


0.5869 - 0.09241 


0.5869 - 0.09241 












(<0.01% , <0.01%) 




1 




0.5715 - 0.28061 


0.5724 - 0.28041 


0.5730 - 0.27981 












(0.10% , -0.21%) 




2 




0.5463 - 0.47511 


0.5462 - 0.47681 




3 







0.7842 - 0.09241 


0.7844 - 0.09241 


0.7844 - 0.09241 












(<0.01% , <0.01%) 




1 




0.7728 - 0.27911 


0.7732 - 0.27901 


0.7734 - 0.27901 












(0.03% , <0.01%) 




2 




0.7524 - 0.47001 


0.7522 - 0.47071 


0.7493 - 0.47751 












(-0.39% , 1.44%) 




3 




0.7257 - 0.66521 


0.7237 - 0.67001 
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TABLE V: The QNM Frequencies for the Spin-\ Perturbation of the Kerr BH with a rotation parameter a = 0.80. 



1 


n 


Numerical 


Third Order WKB(J) 


Sixth Order WKB(J) 


AIM 








0.1932 - 0.0891i 


0.1883 - 0.0896i 

(-2.54% , 0.56%) 


0.1914 - 0.0865i 

(-0.93% , -2.92%) 


0.1920 - 0.0872i 

(-0.62% , -2.13%) 


1 



1 


0.3993 - 0.0893i 
0.3789 - 0.2728i 


0.3956 - 0.0881i 
(-0.93% , -1.34%) 

0.3751 - 0.2701i 
(-1.00% , -0.99%) 


0.3967 - 0.0880i 
(-0.65% , -1.46%) 

0.3777 - 0.2687i 
(-0.32% , -1.50%) 


0.3965 - 0.0880i 
(-0.70% , -1.46%) 


2 



1 
2 




0.5984 - 0.0881i 
0.5844 - 0.2669i 
0.5600 - 0.4512i 


0.5987 - 0.0881i 
0.5855 - 0.2667i 
0.5609 - 0.4517i 


0.5987 - 0.0882i 
(<0.01% , 0.11%) 

0.5847 - 0. 26441 
(-0.14% , -0.86%) 


3 



1 
2 
3 




0.8000 - 0.0882i 
0.7895 - 0.2659i 
0.7702 - 0.4471i 
0.7443 - 0.6320i 


0.8001 - 0.0882i 
0.7900 - 0.2658i 
0.7706 - 0.4473i 
0.7436 - 0.6345i 


0.8001 - 0.08821 
(<0.01% , <0.01%) 

0.7900 - 0.26551 
(<0.01% , -0.11%) 



